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Abstract 

In this brief, an algorithm for controlling chaotic systems using small, continuous time per- 
turbations is presented. Stabilisation is achieved by self controlling feedback using low order LTI 
filters. The algorithm alleviates the need of complex calculations or costly delay elements, and 
can be implemented in a wide variety of systems using simple circuit clmcnts only. 

I. INTRODUCTION 

There has been some increasing interest in recent years in the study of controlling chaotic 
nonlinear systems. The possibilty of obtaining periodic waveforms from a chaotic system 
by stabilising any of the numerous embedded Unstable Periodic Orbits (UPO's) has been 
the guiding control philosophy. The breakthrough in this direction is the OGY algorithm 
[1], which stabilises the UPO by applying ocassional, small, well calculated perturbations 
to the system parameters. A variation of this method is the Ocassional Proportional 
Feedback (OPF) [4], which gives ocassional perturbation to the system parameters, pro- 
portional to the deviation of the system away from the UPO. 

The above methods are in essence discrete in nature, and calculates the perturbation 
based upon local behavior of the system in a neighborhood of the UPO. This is an ad- 
vantage, as knowledge about the entire global dynamics is not necessary. However since 
control is exercised only ocassionaly in a small neighborhood of the UPO, the system 
becomes susceptible to ocassional bursts away from the UPO under moderate noise. To 
overcome this limitations the idea of time continous control was proposed by Pyragas [3] . 
Among the two algorithms suggested by Pyragas one required an external periodic signal 
approximating the UPO, while the other recovered the periodic signal by Delay Coordi- 
nate method. Though various implementations of the second method has been reported, 
the inspiration behind the present report is to replace the delay element which is difficult 
and costly to obtain at some time scales. 

The present algorithm also achieves control by perturbing the system parameters in 
proportion to the error signal between the output signal to be controlled and a periodic 
signal. But the periodic signal is derived from the chaotic attractor itself, by passing the 
chaotic output of the system through a band pass filter with a narrow pass band. The filter 
may be a simple LTI one, which can be implemented using resistor, capacitor and opamps 
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only. It has been observed that though LTI filters are unable to filter out a periodic signal 
from a chaotic one, when connected in the above configuration, the whole system can be 
controlled to a periodic orbit. 

II. CONTROL ALGORITHM 

Let us consider a dynamical system given by the model: 

X = f(X,p) 

(1) 

y = c P x 

Where, p is a parameter available for perturbation and y is the output signal to be con- 
trolled. We assume both y and p to be scalars. The system is connected in the configuration 
shown in Fig 1. 

The filter in the feedback path is band pass type with a narrow pass band. A second 
order notch filter is found to be sufficient for effective control. To select the pass frequency 
of the filter, a FFT of the output signal is obtained, which would be typically spread 
spectrum, with a broad peak centred around frequencies corresponding to the UPO's. We 
select the pass frequency of filter to be within this window. A sufficiently high Q value 
for the filter is selected and assumed to be tunable. If the output of the filter is yj(t) the 
perturbation applied to the system parameter p is of the form 

Ap(t) = K(y f (t) - y(t)) (2) 

The value of gain K is tuned to obtain stabilisation. 

When stabilisation is achieved the output of the filter and the system both become 
periodic and close to each other; and, the perturbation Ap(t) becomes extremely small. 
Therefore, as well as in OGY and Pyragas's method small external force is used for sta- 
bilsation. Also since the pass frequency of the filter was chosen to be in the same window 
as that of the UPO's - the stabilised orbit lies in a small neighborhood of the UPO of the 
original system. The smallness of the perturbation signal depends on how close to UPO 
the system is stabilised and the noise level present. 

Though general validity and sufficiency conditions for the effectiveness of the above 
algorithm is difficult to prove, it appears that stabilisation in the above method is achieved 
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through additional degrees of freedom introduced into the system by the filter in the 
feedback loop. The filter does not change much the projection of the system dynamics 
into the original low dimensional state space, but change only the Lyapunov exponent of 
the UPO. As compared to the delayed feedback method suggested by Pyragas [3], where 
the system dimension increases to infinity, a finite increase in system dimension occurs 
here, i.e. the finite dimensional LTI filter aproximates the delay element in some sense. 
If the state space realisation of the filter is 



(3) 



X f = A f X f + B f u 
Vf = C f X f 

and the system equations for the chaotic dynamics are given by Eqn 1. Then the state 
equations for the augmented system becomes 



X 









(4) 



f(X,p + K(C f X f -C p X)) 
B f C p X + A f X f 

Stability of the UPO's can be studied by linearising the augmented system about a fixed 
point and studying the singular values of the resulting Jacobian matrix. 

A. Stabilisation to higher periodic orbits 

It has been observed that the system could be stabilised to period-2 or higher periodic 
orbits if the frequency components of these orbits are close enough. Stabilisation can be 
achieved in this case by shifting the pass frquency of the filter towards the edge of the pass 
window of the original system and relaxing the Q value. An alternative approach could 
be using parallel combination of filters. 

III. RESULTS FOR SOME GENERIC SYSTEMS 

The Lorentz system is used to illustrate the main results for the algorithm. The system 
equations for the Lorentz system are given by 



Xl = —x\ + x 2 x 3 

x 2 = 3(x 3 - x 2 ) 

X3 = -x±x 2 + rx 2 - x 3 



(5) 
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The state x-i is selected as the output signal to be controlled and r is the system param- 
eter available for perturbation. The nominal value of r is taken as 26.0 . 
The second order filter with transfer function 



F(s) = CU ' 



s 2 + ^s + Ujl 



is used in the feedback loop. The FFT spectrum of the output signal is shown in Fig 2b. 
Accordingly the pass frequency of the filter is choosen as 1.2Hz, which lies in the pass 
window of the chaotic system. Stabilisation of the period-1 orbit is obtained at Q — 8.0, 
K = 0.52 and C = 0.90. Fig 3b. and 3a. shows the stabilsed output signal and the 
perturbation in the system parameter respectively. Fig 3c. shows the FFT spectrum of 
the controlled output. 

It is to be noted that the stabilised period-1 orbit is almost and not truly sinusoidal. 
A small harmonic component is present in the output, though topologically this can be 
considered an period-1 orbit. This leads to the question, whether the objective of control 
of chaos should be stabilisation to sinusoidal outputs, or not sinusoidal but still period-1 
orbits. The answer depends on the application involved. The control algorithm suggested 
above offers a direct control over the harmonic content of the output signal and can modify 
it according to application demand. 

To illustrate the validity of the algorithm for other systems the results for stabilisation 
of the Rossler system to period-1 orbit is shown in Fig 4. 

Local stability of the system can be studied by linearising the augmented state equations 
[4] about the stabilised UPO, and considering the singular values of the resulting Jacobian 
matrix. For the Lorentz system it is observed that the singularvalues become negative for 
K e [0.36,0.55] for Q = 8 and f = 1.2Hz. The variation of singularvalues with K is 
shown in Fig 6b. 

IV. CONCLUSION 

The above algorithm offers in a naive form a paradigm in which the problem of control 
of chaotic systems can also be viewed as that of synchronisation of two back to back, 
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mutually coupled system, of which one may be a chaotic system with numerous embedded 
UPO's and the other may be a LTI or a simpler nonlinear system whose natural dynamics 
is periodic. By synchronisation of outputs of both the systems (by connecting them in 
feedback with a coupling gain) and tuning the simpler system such that its dynamics is 
close to an UPO of the chaotic system, the chaotic system can be stabilised to a periodic 
orbit. It is also observed that the above algorithm is more effective in controlling systems 
having homoclinic attractors, which points towards its similarity with non-linear PLL's. 
The condition on the LTI system for asymptotic synchronisation and stabilisation to UPO, 
is under investigation and will be reported elsewhere. There is scope for future studies in 
using higher order filters and nonlinear systems in the feedback path. 

In conclusion, a method of stabilising chaotic systems approximately to a UPO, by small 
continous time perturbations, is presented. The main advantage of the above method lies 
in its easy applicabilty to a wide variety of systems. No complex calculations are involved 
and the algorithm can be implemented using low order LTI filters only. 
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Fig. 1 

Block Diagram of the Control Algorithm 
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Fig. 2 

Dynamics of the uncontrolled Lorentz system, (a) Phase potrait and output (x 2 ) for 

THE UNCONTROLLED SYSTEM, (b) FFT SPECTRUM OF THE OUTPUT OVER A TIME SCALE OF 50s, 
SAMPLED AT 50Hz, (c) FREQUENCY RESPONSE OF THE FILTER USED IN FEEDBACK PATH. 
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Fig. 3 

Controlled period- 1 orbit for the Lorentz system: (a) Phase Potrait and perturbation 
in r, (note that unlike Pyragas's method where perturbation can be large during the 

TRANSIENT PERIOD, LEADING TO PROBLEMS LIKE MULTISTABILTY, HERE PERTURBATION IS SMALL 
EVEN DURING THE TRANSIENTS), (b) OUTPUT SIGNAL, (c) FFT SPECTRUM OF THE OUTPUT AFTER 
THE TRANSIENTS ARE OVER, (THE PEAK AT DC IS DUE TO THE DC BIAS OF THE OUTPUT), THE 
STABILISED ORBIT IS AN almost P-l ORBIT, A SMALL HIGHER FREQUENCY COMPONENT IS PRESENT 
ALONG WITH THE MAIN FREQUENCY COMPONENT AT / « 1.2Hz. 
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(b) 



Fig. 4 

Controlled dynamics for the Rossler system: 
x\ = -x 2 - x 3 ,x 2 = x\ + 0.2x 2 + Ap(t), x 3 = 0.2 + x 3 (xi - 5.7), perturbation Ap(t) is determined 
by Eqn.2. with y(t) = x 2 (t) (a) Time plot of Ap(t) (b) Controlled output, (c) FFT 

SPECTRUM OF THE UNCONTROLLED OUTPUT, (d) FFT SPECTRUM OF THE CONTROLLED OUTPUT. 

Control is achieved at f = 0.12Hz, Q = 1.8, C = 1.0, K = 0.473. 
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Fig. 5 

(a) Poincare section for the uncontrolled Lorentz system OBTAINED AT X\ = 25.0 , '+' 
shows the poincare section for stabilised period- 1 orbit, note that the stabilised point 
lies close to the fixed point of the original system (which appears denser in the 
figure), (b) Singular values of the Jacobian matrix for the augmented system taken 
about the period- 1 fixed point shown in poincare section. plots are drawn for the 
variation of singular values over k, for / = 1.2hz and for q = 6,8,10. respectively. 
The value of K is to be chosen in the region where the singularvalues are negative. 
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